Symmetry breaking and restoring wave transmission in diode-antidiode double chains 
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We introduce a system of two parallel-coupled discrete nonlinear Schrodinger (DNLS) inhomoge- 
neous chains. Each one favors the unidirectional transmission of incident packets, in the opposite 
directions with respect to each other. Two different configurations of the diode-antidiode pair 
are considered, a ladder and a plaquette. They feature, respectively, the uniform transverse lin- 
ear coupling, or the coupling limited to the central nonlinear segment of the system. In the case 
of weak linear coupling, the symmetry breaking is observed (i.e., the pair still features the diode 
behavior), while the moderately strong coupling restores the symmetry, making the transmission 
effectively bidirectional. In the case of the ladder, an oscillatory dependence of the transmission on 
the strength of the coupling is observed and qualitatively explained. 

PACS numbers: 05.45.-a 63.20.Ry 



I. INTRODUCTION 

Asymmetric (nomeciprocal) wave propagation induced 
by nonlinearity is of great interest as a basic principle for 
the energy- flow control. It emerges in several different 
context. Among the first examples discussed in the lit- 
erature is the asymmetric phonon transmission through 
a nonlinear interface layer between two very dissimilar 
crystals [1]. In the realm of acoustics the possibility of re- 
alizing a wave diode has been demonstrated for nonlinear 
phononic media [2, 3]. Another promising context is the 
propagation of acoustic pulses through granular systems. 
In the latter context, experimental studies have demon- 
strated a change of the reflectivity of solitary waves from 
the interface between two granular media Q. More re- 
cently, demonstration of rectification of the mechanical- 
energy transfer at sonic frequencies in a one-dimensional 
array of particles has been reported too [E'] . In the field of 
nonlinear optics, the so-called all-optical diode was first 
proposed in Ref. [gI, , and later on realized experimen- 
tally [H. Other diverse implementations of the unidirec- 
tional transmission have been suggested, ranging from 
left-handed materials [9|, periodic [10] and quasiperiodic 
[m photonic crystals, coupled nonlinear cavities [l2| to 
T^T— symmetric waveguides [l3|, An extension to 

quantum regimes, in which few-photon states may also 
display a diode effect, has been proposed in Ref. [15|- A 
related work for electric transmission lines was reported 
in Ref. [iG^. Interestingly, for nonlinear lattices brought 
in contact with thermal baths reciprocity violations are 
related to the second principle of thermodynamics p^Tf . 
Recently, a critical discussion of the topic, including the 
case of thermal sources, was given in Ref. [18]. 

A reference model to discuss this type of phenomena is 
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the Discrete Nonlinear Schrodinger (DNLS) equation [19| 
with embedded nonlinear elements |2Q|-[2^. Beyond its 
relevance in various physical contexts, ranging from non- 
linear optics to the physics of trapped cold gases [25|, [26| , 
the DNLS equation has the major advantage of being one 
of simplest dynamical systems amenable to a complete 
theoretical analysis. In the present context, DNLS-based 
models are particularly relevant, as they allow one to 
solve the scattering problem exactly, without complica- 
tions of dealing with wave harmonics [l^. In particular, 
it has been demonstrated in Ref. that the effect of 
the nonlinear shift of a resonance can be exploited to de- 
sign a diode chain, capable to transmit waves with the 
same amplitude and frequency differently in opposite di- 
rections. 

The propagation direction admitted by a diode chain is 
reversed in a chain which is a specular counterpart of the 
given one. This obvious fact suggests us to consider the 
transmission of waves in double chains, formed by two 
parallel diode cores with opposite orientations, which are 
linearly coupled to each other in the transverse direction. 
This setting, which may be called a diode-antidiode pair 
(ahas a diode dipole), is the subject of the present work. 
The transmission of waves in either direction in such a 
system obviously implies a spontaneous breaking of the 
antisymmetry between the transversely coupled counter- 
oriented cores. It is relevant to mention that the sponta- 
neous symmetry breaking in a system of linearly coupled 
identical one-dimensional lattices with the onsite cubic 
nonlinearity, which admit the transmission of waves in 
both directions, was studied in Ref. [28|. That system 
was modeled by two linearly coupled DNLS equations, 
and the symmetry breaking was considered for quies- 
cent discrete solitons. The spontaneous (anti) symmetry- 
breaking effect, addressed in the present work, is com- 
pletely different, being dynamical, rather than static, and 
realized in terms of incident, transmitted, and reflected 
linear waves, rather than solitons. In this respect, this 
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is a new setting for the study of spontaneous-symmetry- 
breaking effects, which manifest themselves in a great 
variety of nonhnear systems [29j . 

The rest of the paper is organized as fohows. Two 
different models of the diode-antidiode pair are formu- 
lated in Section ini In Section UTTl we describe simulations 
of the wave-packet scattering. In Section |lVl we reports 
results for various dynamical regimes, focusing on the 
spontaneous symmetry breaking and restoration for the 
two different types of the coupled diode-antidiode sys- 
tem. We complete the paper with concluding remarks in 
Section IVl 



We assume the usual scattering setup which implies 
that Un-, ocn and j3n in Eqs. ([1]) are nonvanishing 
only in a finite segment of the chain. 



l<n<N. 



(4) 



The two semi-infinite portions of the lattice (n < 1, 
n > N) represent two uniform linear leads, where waves 
propagate freely. In the ladder configuration, the free 
propagation obeys the following dispersion relation for 
wavenumber k and frequency uj: 



uj{k) 



-2 cos k ± k| 



(5) 



II. MODELS 

We consider two linearly-coupled DNLS chains, with 
onsite dynamical complex variables. Unit) and Vn{t)^ 
where n is the discrete coordinate in the one-dimensional 
chains, and the overdot stands for the time derivative: 

iiln = UnUn-Un^l-Un-l^an\Un\'^Un^ l^nVn, 
iVn VnVn - Vn^l - Vn-1 ^ f3n\Vnfvn ^ l^nUn, (1) 

The chains are non-uniform, with the real onsite non- 
linearity coefficients, an and /3n, as well as the real 
transverse-coupling coefficient, /^n, depending on n. In 
addition, each chains carries a discrete linear potential, 
represented by real discrete functions Un and Vn. 

As sketched in Fig. [TJ we will consider two versions of 
the general model. The first system includes the uniform 
linear coupling, while in the other one the coupling acts 
only at two central sites (with n = 1 and n = 2): 



Kn = for every n. 



(2) 
(3) 



We will refer to these two models as the ladder or pla- 
quette , respectively. In fact, the ladder configuration for 
static discrete solitons in the usual bidirectional system 
was introduced in Ref. [iH, while the plaquette resem- 
bles the bidirectional system with the coupling applied 
at the single central site (soj . 
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the respective group velocity being 

Vgr = duj/dk = 2 sink. 

Actually, we consider the dimer, with TV = 
O, i.e., 

Un = Vo[{l^e)5n,i^{l-e)5n,2] 

Vn = Vo[{l-€)6n,l^{l^e)6n,2]^ 
an= Pn= ^n,l + ^n^- 



(6) 

2 in Eq. 

(7) 
(8) 
(9) 



The system is invariant with respect to the mirror reflec- 
tion about the dimer's center, n ^ A/" + 1 — n, combined 
with interchange of the chains {un^Vn) (vn^Un)- 

The dimer in the single chain {hc = 0) was considered 
in Ref. Numerical simulations have demonstrated 

that Gaussian wave packets impinging on the dimer from 
the two opposite directions have very different transmis- 
sion coefficients, which implies that the single chain with 
the embedded dimer plays the role of the diode for the 
wavepacket transmission. 



III. SIMULATIONS OF THE WAVEPACKET 
SCATTERING 



We performed simulations of Eqs. ([T]) for the chains 
of finite lengths, i.e. for \n\ < M (this means that each 
of the two chains is composed of 2M + 1 sites). Open 
boundary conditions are enforced on both chains. Initial 
condition were taken as two Gaussian wavepackets, with 
no < 0: 



Un{0) = /(l + (^)exp 
Vn{0) = ±I{1-S)exp 



(n - no)^ 



(n - no)^ 



- ikon 



ikon 



,(10) 



FIG. 1: (Color online) A sketch of the two models (for = 2): 
the "ladder" in the upper panel, and the "plaquette" in the 
lower one. In both models, the nonlinearity acts at two central 
sites of the chains. 



where ^ is a small perturbation (typically we set ^ = 10~^ 
henceforth), which seeds symmetry-breaking effects. The 
symmetries of the models suggest to consider two types 
(even and odd) of initial conditions (p!Q|) . with plus and 
minus signs, respectively, in the expression for Vn{0). 
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Note that the total initial input power [which is a dynam- 
ical invariant of Eqs. O] is (|iin(0)P + |vn(0)P) (X 
w\I\'^ (up to a small correction ~ S). 

The wavepacket-transmission coefficients in the two 
coupled chains, produced by the simulations in the course 
of temporal interval < t < tfin, are naturally defined as 

" E„<i(K(o)P + K(o)|2)' ^ > 
""E„<i(K(0)|2 + K(0)|2)- ^ > 

To minimize the dispersive effects and, thus, the depen- 
dence of the scattering process on the initial position no, 
we fix /co = 7r/2 in Eqs. (p!Q|) . Moreover, for a given lat- 
tice size M, tfin is chosen to be not too large, so as to 
avoid hitting the boundary sites by the transmitted and 
reflected packets (M = 500 and tfin = 250 are the values 
used in the following). 

First, Fig. [2^ shows that, for the case of very weakly 
coupled chains, with = 0.01, the transmission of the in- 
cident packet is strongly different in the two components 
{tu = 0.34, ty = 0.15), as expected from the analysis of 
the single chain [g^ (in that work, the parameters were 
the same as used there). However, even a moderate cou- 
pling, = 0.1 (Fig. [2]3), restores the symmetry between 
the chains, giving tu ^ ty. 

Figure [3] shows a snapshot of the norm densities at the 
end of the same run which is presented in Fig. [JJd. Apart 
from the distortion of the reflected and transmitted pack- 
ets, it can been seen that some energy remains trapped 
by the central segment in both cases, being larger in the 
first one (see insets). A similar phenomenology was ob- 
served in the single chain [3l| , which is explained by the 
spontaneous formation of localized defect modes at the 
nonlinear sites (see again the insets in Fig. [3]). The scat- 
tered packets remain almost monochromatic, at the in- 
cident wave number /cq, with some weak radiation waves 
leaking throughout the lattice. 
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FIG. 2: (Color online) Numerical simulations of the trans- 
mission of Gaussian wavepackets ([T0|) in the ladder system. 
Parameters are Vo -2.5, A;o = 7r/2, e = 0.05, M 500, 
|/|^ = 3, w = 100 and no = —250. Panels (a) and (b) pertain, 
respectively, to the weak and moderate couplings, k — 0.01 
and HI — 0.1. 



A. The ladder-type system 



IV. SYMMETRY BREAKING AND 
RESTORATION 

The results of the previous Section indicate that the 
symmetry of the scattered packets depends on the cou- 
pling. In the present Section we present a more detailed 
analysis of the transition from symmetry-broken to sym- 
metry restored outputs as a function of the parameters. 
We anticipate that the scenario may be considerably 
complex, as the transitions from the two regimes may 
sensitively depend not only on the interchain coupling 
but also on the nonlinearity level. 



In Fig. m we show the transmission coefficients, de- 
fined according to Eqs. (pT]) and (fT2]) . as functions of the 
wavepacket's intensity, For the very weak coupling, 
there is a region in which the two transmission coeffi- 
cients are significantly different (as in the example shown 
above), i.e., the diode-antidiode pair still behaves similar 
to a single diode chain. However, the difference between 
tu and ty tends to disappear with increasing n. For in- 
stance, at = 1 the outputs are almost indistinguishable, 
even at the central sites, which implies that the suffi- 
ciently strong coupling makes the transmission effectively 
isotropic (bidirectional). The results strongly depend on 
width w of the incident packet, see Eq. (p!Q|) . For small 
w the curves are much smoother, while broader pack- 
ets (e.g., with w = 100) give rise to wilder oscillations. 
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FIG. 3: (Color online) The snapshot of the norm densities in 
the two chains at the end of the simulations, for the ladder 
with moderate coupling, k = 0.1. The insets zoom the fields 
around the nonlinear sites, showing that some energy remains 
trapped after the scattering. 



Qualitatively speaking, this can be understood by noting 
that the incoming wave packet acts as a time-dependent 
force on the nonlinear core. Increasing w amount to an 
increase the power and thus the nonhnearity, resulting in 
more complex (possibly chaotic) dynamics, which yields 
large changes in the transmission and thus a sensitive 
dependence on the parameters of the incident packet. 

A noteworthy phenomenon is shown in Fig. [5l where 
we display the transmission coefficients as a function of 
the interchain coupling for two fixed intensities, 
Namely, an antiphase oscillating behavior of tu and ty is 
observed for the smaller value of 

To understand the origin of such oscillations, in Fig. 
[6] we report the limit case of the linear lattice, with 
<^n = /^n = in Eq. ([T]). In this case, the same os- 
cillations occur too, although with a small amplitude. 
Thus, the oscillations of the transmission coefficients in 
the linear ladder are enhanced by a moderately strong 
nonlinearity. An essential point is that this effect occurs 
in the ladder but not in the plaquette. This fact, as well 
as the entire effect of the oscillations, is explained by the 
fact that, in the ladder, the waves in the linearly coupled 
top and bottom chains may periodically switch between 
them after the passage of the core segment (the two cen- 
tral sites). 

To analyze the oscillations in a quantitative form, we 
note that, according to Eqs. ([T]), the temporal frequency 
of the energy exchange between the chains is In the 
combination with the group velocity given by Eq. ([6]), 
this implies that the spatial period of the exchange in 
the wave packet propagating in the uniform linear lad- 
der is L = 27r]^r//^ = 4:7rK.~^ sink. Thus, the intensity 
of the transmitted wave oscillates between the chains. 
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FIG. 4: (Color online) The dependence of the wavepacket 
transmission coefficients tu (solid lines) and tv (dashed), see 
Eqs. (|lip and ()12p . on the intensity |/|^ for the ladder model 
at two different interchain couplings Other parameters are 
as in Fig. O 




FIG. 5: (Color online) The dependence of the wavepacket 
transmission coefficients tu (solid lines) and tv (dashed), see 
Eqs. (pT]) and (p^ . on coupling k for the ladder model at two 
different intensities. Other parameters are as in Fig. [2] 



as a function of discrete coordinate n as cos(27rn/L) = 
cos ^(2sinA:)~^ nn^. At the edge of the lattice, n = M, 
this implies oscillations between tu and t^, following the 
variation of z^, in the form of cos ^(2sin/c)~^ kM^. The 
period of the latter oscillations is 



Ahi 47tM~^ sink. 



(13) 



For the value k = 7r/2 fixed throughout the paper and 
for M = 500 (see caption to Fig. Ej), Eq. 1^ yields 
Ak, = 47r/500 ~ 0.025, which agrees with the picture 
observed in Fig. [6l 

To identify the symmetry-broken and symmetry- 
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FIG. 6: (Color online) The dependence of the wavepacket 
transmission coefficients tu (solid lines) and tv (dashed), see 
Eqs. and (|12p . on coupling k for |/|^ = 0.5 in the linear 
system, with an = I3n = 0. 



restored regions, we define the asymmetry factor, 



til, ~\~ 



(14) 



which is small when the output intensities in the two 
chains are nearly equal, and approaches unity when one 
of the two transmission coefficients vanishes. 

In Figs. and [8] we show contour plots of a in the 
(|/p,Ai:) plane for different wavepacket widths w in the 
ladder model with even and odd initial conditions, re- 
spectively, in Eqs. (fTQ|) . It is seen that there are pa- 
rameter ranges where a is very small, alternating with 
many regions where the symmetry breaking is sizeable. 
Roughly, the regions of small a are concentrated around 
|/p ^ 0.6 (in the leftmost part of the figures) and 
for larger intensities and couplings — around |/p > 2.5, 
n > 0.02. Notably, there are oscillations caused by the 
antiphase behavior observed in Fig O 



B. The plaquette-type system 

As seen in Fig. [9l the structure of symmetry-broken 
regions is much simpler for the plaquette configuration. 
This is explained by the fact that the absence of the 
linear coupling between the chains outside of the central 
segment does not allow more subtle feature to develop 
via the continuing switching of the power between the 
chains. 

In the case of the plaquette, it is also possible to launch 
an input in a single chain, as the chains are decoupled 
outside of the central segment. In this vein, we have 
investigated the asymmetric initial conditions with the 
wavepacket initially launched in one of the two chains: 



^n(O) = /exp 

^n(O) = 0, 



(n - no)^ 
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FIG. 7: (Color online) The ladder configuration with even 
initial conditions, given by Eqs. ([T0|) with plus sign. Contour 
plots of the asymmetry factor, (p^ . are shown for different 
wavepacket's widths, w = 10, 20, 30, 100, from top to bottom. 
Other parameters are as in Fig. [2l 




(15) 
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FIG. 8: (Color online) The ladder configuration with odd 
initial conditions, given by Eqs. (p^ with the minus sign. 
Contour plots of the asymmetry factor, ([T4j), are shown for 
different wavepacket's widths, = 10, 20, 30, 100, from top to 
bottom. Other parameters are as in Fig. [2] 
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In the case of the ladder, an oscihatory dependence of the 
0.08 



0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
iii2 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
iii2 



FIG. 9: (Color online) The plaquette configuration, with the 
even initial conditions, given by Eqs. (p^ with the plus sign. 
Contour plots of the asymmetry factor, ([H]), are shown for 
different wavepacket's widths, w = 10, 20, 30, 100 from top to 
bottom. Other parameters are as in Fig. [2] 



as well as 



Un{0) = 0, 

Vn{0) /exp 



(n - no)^ 



- ikon 



(16) 



For comparison with the previous cases, we consider 
the same ranges of coupling parameters and intensities. 
Here, the results are rather simple, as shown in Fig. [TOl 
The growth of the intensity, /, i.e., of the nonlinearity 
strength, causes the wave to stay in the original chain, 
while the increase of the coupling constant, z^, pushes the 
system towards even distribution of the intensity between 
the coupled chains. 



V. CONCLUSIONS 

We have presented numerical simulations for the scat- 
tering of wavepackets in a model of two parallel-coupled 
lattices, each of them acting individually as diode chains 
in opposite directions. We discussed two different con- 
figurations, the ladder and the plaquette. In the case of 
weak linear coupling, the symmetry breaking is observed 
(i.e., the pair still features the diode behavior), while the 
moderately strong coupling gradually restores the sym- 
metry (makes the transmission effectively bidirectional). 



FIG. 10: (Color online) The plaquette configuration with ini- 
tial conditions Eq. (|15p corresponding to the wavepacket ini- 
tially located in the single chain. Contour plots of the asym- 
metry factor, (O, are displayed for different wavepacket's 
widths, w = 10, 20, 30, 100 from top to bottom. Other pa- 
rameters as in Fig. [2l 



transmission on the strength of the coupling is observed. 
Such oscillations are qualitatively explained as nonlin- 
ear versions of those occurring in the linear limit. As 
a result, the symmetry-broken regions in the parameter 
space form a complex pattern, alternating with regions 
of the restored symmetry. The location of such regions 
are also strongly dependent on the wavepacket widths. 

As mentioned in the Introduction, one our motivations 
was to study asymmetric wavepacket propagation p7| . 
In this context, the present study is a first attempt to 
discuss a "multichannel" configuration in which several 
inputs can be simultaneously managed through nonlinear 
interactions. We outlined how the concepts of symmetry 
breaking and restoration are important issues. Indeed, 
they directly relate to the functionality of the proposed 
system as a tool to selectively control the transmission 
direction of multiple inputs. Moreover, the present pa- 
per is the a starting point to address reciprocity viola- 
tions in general lattices with nontrivial topologies and/or 
couphngs. This work can be extended by considering 
more complex setups where multiple "leads" are coupled 
through nonlinear sites. In particular, it may be inter- 
esting to consider a network built of alternating parallel- 
coupled chains with opposite transmission directions. 
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